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Abstract 

Motivated by the precision results in the electroweak theory studies of two-loop Feynman diagrams 
are performed. Specifically this paper gives a contribution to the knowledge of massive two-loop self- 
energy diagrams in arbitrary and especially four dimensions. This is done in three respects: firstly 
results in terms of generalized, multivariable hypergeometric functions are presented giving explicit 
series for small and large momenta. Secondly the imaginary parts of these integrals are expressed as 
complete elliptic integrals. Finally one-dimensional integral representations with elementary functions 
are derived. They are very well suited for the numerical evaluations. 
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1 Introduction 



The beautiful results of the LEP1 experiments have shown that the electroweak theory has a predictive 
power like that of QED several decades ago. It is to be expected that eventually the electroweak theory 
will provide high precision predictions for many experiments in the present and near future. One may in 
particular think in this respect of the measurement of the Z mass and width, a better determination of 
the W mass and a first indication of the value of the top mass. This will in the future require two-loop 
calculations in the electroweak theory. One then has to face several problems: amongst others, a large 
number of diagrams and many different non-negligible masses. It is known that analytical results for arbi- 
trary non-vanishing masses are in general not obtainable for two-loop diagrams in the form of generalized 
polylogarithms. The present paper contributes to the knowledge about this problem by introducing and 
applying new analytical and numerical approaches. A number of new results will be presented. 

We focus on the simplest class of two-loop diagrams, i.e. the set of scalar two-loop self-energies. They 
are essential for physical predictions and moreover they show the typical problems one encounters in the 
evaluation of two-loop diagrams. There exist four non-trivial self-energy diagrams 
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Figure 1: Topologies 



In order to classify our conventions we give the explicit form of the so-called master diagram: 

1 1 



Ti2345(p ;m 1 ,m 2 ,m 3 ,m 4 ,m 5 ) =« 



» 



mf 



where each bracket denotes a D dimensional integration over one loop momentum q: 

d D q 



<...>= 



iTT 2 (2lTfl) 



D-4 



(•••), 



(1) 



(2) 



where fi is an arbitrary mass. The momenta of the propagators ki,...,ks are determined by momentum 
conservation in terms of the external momentum p and the integration momenta qi and q 2 . It can be 
immediately seen that the diagram Tii234(p 2 ; m\, mf, m 2 ,, m§, m|) can be reduced to a difference of two 
T1234 diagrams by partial fractioning of the propagators with the same momentum: 



mf 



,/2 



(T 1234 (p 2 ;mlmlmlml) 



(3) 
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In the special case mi = m! x this has to be taken as a derivative with respect to m 2 . So we are only left 
with three diagrams to be discussed. Of these the master diagram is convergent in four dimensions, the 
other two are divergent. After evaluating them in D = 4 — 25 dimensions, one would like to derive an 
expansion in S which starts in 1/S 2 and 1/8 terms. 

In recent years several groups derived small and large p 2 expansions for these diagrams. For the 
general massive case they obtained algorithms for symbolic evaluation of these expansion coefficients. For 
the general coefficients in these expansions, no closed expression was derived. More recently multiple 
series were found for which the general coefficient is known and which lead to generalized hypergeomctric 
functions in several variables. One aim of this paper is to extend these results, an other to express the 
imaginary parts by elliptic integrals. On the other hand an elegant numerical method has been worked 
out leading to a two-dimensional integral representation 0, ^ . The second aim of our paper is to improve 
on this by introducing a one-dimensional integral representation containing only elementary functions. 
For this a general procedure for diagrams containing one-loop self-energy expressions in a subloop is 
presented. The actual outline of the paper is as follows: in sect. 2 the analytic approach leading to 
generalized hypergeometric functions and the corresponding explicit series representations are given, the 
elliptic integrals for the diagrams are calculated in sect. 3. The next sect, describes the numerical 
approach by means of a one-dimensional representation which is derived by a self-energy insertion into 
one-loop diagrams. Finally sect. 5 presents some numerical comparisons and draws conclusions. 

2 Analytic approaches and hypergeometric functions 

In this section we present analytic results in an arbitrary number of dimensions D for the two-loop scalar 
self-energy diagrams. They will be in the form of generalized hypergeometric functions, that is multiple 
series of ratios of p 2 and the masses. We also discuss the expansion of these results around D = 4. 

The simplest case of a two-loop scalar self-energy diagram Ti23(p 2 ;m 2 ,m 2 ,m 2 ) was discussed in a 
previous paper In fact, the result turned out to be represented by Lauricella functions which could 
be derived by x-space techniques or Mellin-Barnes representation for a massive propagator. It turns out 
that there is a third method using dispersion relations. Since the latter method will be used several 
times in this paper we shall first rederive the result for T 12 3(p 2 ; raf, m|, m§) by means of dispersion re- 
lations. Then we focus our attention on the two-loop scalar self-energy diagram with four propagators 
Ti234{p 2 ', m 2 , m|, m§, m 2 ) where the small p 2 behavior can again be derived by the dispersion method and 
the large p 2 behavior by the Mellin-Barnes representation method. 

Although we have three methods, the master diagram still poses many problems, the origin of which 
we will briefly discuss in the last section. 

So we start with the rederivation of T\23(p 2 ; m 2 , m 2 , m 2 ) by means of dispersion relations. The imagi- 
nary part is given by 
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(\/p 2 -mi) 



ds 



A 2 d (s,W2,m|)A2 6 (p 2 ,s,mf) 



,1-6 



(m 2 +m 3 ) 2 



2M-6 



(m2+m.3) 2 

where A(a, 6, c) = (a — b — c) 2 — 46c is the Kallen function. The dispersion relation reads 



rrt { 2 2 2 „,2\ 

Ti23(P ;m 1 ,m 2 ,m 3 ) 



1 

(mi+m 2 +m3) 2 



dz- 



rAT^^m 2 ,™ 2 .,?™ 2 ) 



Z — p* 



1 

4^2 



ds ABq(s; m 2 , m\) 



dz 



(m2+m 3 ) 2 



z — p* 



ABo{z;s,ml) 



Using the expansion 



I 1 fp 2 



z — p* z 



V 5^ I y 



k=0 



we perform first the integration over z: 

oo °° 

A = $> 2 ) fe ( dz z^-^z-uf/^iz-v) 1 / 2 - 



k=0 



= J2( p2 ^ u ~ S ~ kB ( k + 5 ^/ 2 ~ 5 hFi{5 - 1/2, k + S; k + 3/2; v/u) . 



k=0 



(4) 



(5) 



(6) 



(7) 



where u = (mi + t/s) 2 an d v = ( m i ~ V 7 ^) 2 - One transforms the Gauss hypergeometric function 2 Fi 
using relations which one can find in Q 

2^(5 - 1/2, k + 6; k + 3/2; v/u)= 



2yfs 



-2(k+S) 



mi + y/S / 

r(i-6) 



+ 



T(k + 2-S) 

m 2 V' A " T(6-l) 
T(k + S) 



2 F 1 (k + 5, k + 1; 2k + 2; (1 - m 2 /s)) 
2 F 1 {k + 5,k+l-8;m 2 1 /s) 



2 F 1 (k + 2-6,k + l;2-S;m 2 1 /s) 



2yfs 



mi 



-2(k+S) 



r(2fc + 2) 

jfe! 



(8) 
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With this result, Ti 23 (p 2 ;m\,m\,m z v ) becomes: 

T\2z{p 2 ;rn\, ml, ml) = 

o. s T(l-S)r(8-l) 



-(47T/T 



T(2 - 28) 



l \--\s,mi,mi) 



(m2+m,3) 2 



*? ( 



TOT / TOT 



-S 



Fi(l,2-S;2-S,2-5; p 2 /s, m\/s) 



F 4 (l,8;2-8,8;p 2 / S ,m 2 1 /s) 



(9) 



Since in (^|) the z integration up to some factors represents the one-loop self-energy, we find as a byproduct 



B (p ;m 1 ,m 2 ) 



T(8-l) 



m\ ( m\ 
m 2 \ 4tt/i 2 



F 4 (l,2-8; 2-8, 2-8;p 2 /m 2 , m 2 /m 2 ) 
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F 4 (M;2- 8,5;p 2 •/mlml/ml) 



{i- Xl )- d {i-x 2 y 



1-8 

+(xi ' a*) 1-2 ^! - 8,1 - S)(-l)- s } 



2 F 1 {l,S,2-8: 



1 - Xi , 
1 — Xl' 



where 



xi,2 = ^2 \ P -m 1 + m 2 ± \l \{p 2 ,m\,m 2 2 ) ) . 



(10) 



(11) 



Using the definition of the F 4 functions we easily perform the integration over s. After some manipu- 
lations similar to the ones used in (^) we obtain the result for the London transport diagram which agrees 
with the one derived using x-space techniques or Mellin-Barnes representation 



Ti2s{p 2 \rn\,ml,ml) = -m\ 



( mj 



2(i/-l) 



\z\zl T 2 {-v) F<i 3) (l, 1 + v; 1 + v, 1 + v, 1 + v; z x ,z 2 , z 3 ) 

-z\ T 2 {-v) F ( c\l, 1 - v; 1 + v, 1 - v, 1 + v; Zl ,z 2 ,z 3 ) 
—z 2 T 2 (-v) Fg\l, 1 — v;l — v,l -\- v,l + v\ z 1 ,z 2 ,z 3 ) 
-Y{y)Y{-v)V(l - 2v) Fg\l - 2v,1 - v;1 - u,1 - u,1 + p; z ly z 2 , z 3 )} 



(12) 
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where z% = mf/m§ , i = 1, 2 ,23 = p 2 /m 2 and 1/ = 1 — 8. The above result is valid for small p 2 but it 
can be analytically continued to the large p 2 region, using known analytic continuation formulae for the 



Lauricella function F. 



(3) 



c 

With the dispersion method described above we derive the small p 2 result for X1234 in D dimensions. 

(2) 

The discontinuity AT 123 4 is a sum of a two and a three-particle cut, which we denote by AT^ 34 and 
AT 12 3 4 , respectively. The two-particle cut is given by 

,(2) 



AT^(p 2 ;m 2 ) = AB (p 2 ; ? t 1 2 , m 2 )B (m 2 ; m 2 , to 2 ) , (13) 

where according to the Cutkosky rules m 2 in Bo(m 2 ; m 2 , m 2 ) is considerd as to 2 + ie, when mi > TO3 + UI4. 
Inserting this result in the dispersion integral gives 



1 



dz ■ 



1 



z — p" 



ATgl 4 (z;m 2 ) 



(mi+m 2 ) 2 

B (p 2 ; to 2 , m 2 ) B (to 2 ; m 2 , to 2 ) 



(14) 



The discontinuity AT 42 3 4 is related to the three particle-cut given in (^) but with an additional factor, 
which according to the Cutkosky rules is the complex conjugate of the propagator (s — to 2 + ze) _1 and 
therefore we get the dispersion relation 



.2. ■■> , , 1 ,,,2\2<5 r 2 (i 6) 



S(p>i) = -(V)' 



T 2 (2- 25) 



d.s 



dz 



(m 3 +m 4 ) 2 

\$- 5 (s,z,m%) 
z 1 - s (z-p 2 ) 



A2 5 (s,m 2 ,ml) 
s 1 ^ s (s — m\ — ie) 



(15) 



(2) (3) (2) (3) 

A further discussion of the analytic properties of AT 12 g 4 , Ar i2 g 4 and consequently T 12 3 4 , T 12 3 4 is postponed 
to subsection 3.2. After performing the z integration we obtain: 



^ 4 ( P 2 ;m?) = -(4 7 TM 2 ) 5r(1 - W - 1) 



T(2-25) 



ds- 



7X 2 (3,7713,7714) 



(ui3+m4) ! 



S \4:TT[1 2 
-5 



F i (l,2-S;2-S,2-5;p 2 /s,m 2 2 /s) 



\47T/J 2 



F 4 (l,<5;2-5,,5;7 J 2 / S ,m 2 / S ) 



(16) 



Expanding (s — to 2 ) 1 in m 2 /s and performing the integration over s one gets 



T(l-6)T(l + 6) f m\ 

5 \47T/J 2 



-28 
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k nn \l m T(l + m + n)r(l + 6 + m + n + k + l) 
x ^ ^1^(1-^3)^4 r(2-<S + 7n)m!nm 

1 _ <5 r(2 - (5 + m + n)r(2 + m + n + fc + Z) 
^ 2 r(2 - 6 + n)T(2(m + n + k)+4 + l) 

T(6 + m + n)T(2S + m + n + k + l) 
' T(S + n)r(2(m + n + k) + 2 + 26 + 1) 



(17) 



where = m 2 /m| with i = 1,2,3 and 24 = p /m\. Note that the contribution from the three-particle 
cut is written in terms of multiple series which are no longer Lauricella functions but rather belong to a 
special class of generalized hypergeometric functions. To our knowledge they have not been studied in the 
mathematical literature. One can however obtain information on the convergence region. The series in 
(1 — Z3) can be written as 2-F1 functions which can be transformed into 2-P1 functions with variable Z3. One 
then obtains four quartic series in Zj. Applying the standard reasoning (see e.g. 0) we get the following 
conditions for convergence 

m 2 + \/\p 2 \ < nil and mi + m,3 < 1114. (18) 
With these results T1234 becomes 

Ti234(p 2 ;m 2 ) = B Q (p 2 ;ml,rnl)B (ml;rnj,ml) + T^l 4 (p 2 ; mj) . (19) 

Next the general small p 2 result for T1234 will be expanded in 6. 

The following combination of the general massive case with a massless case is chosen in such a way 
that the infinite parts cancel B 

Tl234Ar(p 2 ;TOi,Tl2: m 3\ TO 4) = T 12 3i{p 2 \ m\, m\, m\, ml) 

-T 1234 {p 2 ;mj : mlO,0). (20) 

An analytic form of this combination is obtained by expanding the multiple series and their coefficients 
in 5, where the first and the second logarithmic derivatives of the r-function occur at integer arguments 

n n 1 

V»(n + l) = -7 + Ep ^'( n + !) = C( 2 ) -Es- ( 21 ) 

k=l k=l 

with the Euler constant 7 and £(2) = it 2 /6. 

The l/^ 2 , 1/(5 and 7 terms indeed drop out from the result and a finite combination of various mul- 
tiple series remains and some other terms corresponding to the finite parts of Ti234(p 2 ; m 2 , m|, 0, 0) and 
Bq(p 2 ; m 2 , TO2)-Bo(to 2 ; m|, m 2 ) which are given in 0. The result is given in the appendix. 

The iarge p 2 result has been derived using the Mellin-Barnes representation for a massive propagator 

+200 

1 1 1 f (-m 2 ) s 



{k 2 - to 2 )" r(a) 2ni J {k 2 ) a + s 



ds 7T^n-m^ + s), (22) 
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where the integration contour in the s plane must separate the series of poles of r(— s) on the right from 
the series of poles of F(s + a) on the left. In the expression for T1234 we apply ( |2^ ) with a = 1 to all 
propagators, thereby relating the general massive case to the massless one, but with the propagators raised 
to arbitrary powers 1+Sj. The required expression is well-known, see e.g.H 

— in Si 4 

Tl234(p 2 ; 0, 0, 0, 0; 1 + 81, 1 + S 2 , 1 + S 3 , 1 + S 4 ) = e (4 ^ 2) D-4 (-pY^^ " 

r(2 + .s 3 + s 4 - D/2)T(D/2 - 1 - s 2 )T(D/2 - 1 - s 3 )r(D/2 - 1 - s 4 ) 

r(i + S2 )r(i + S3 )r(i + s 3 )t(d - 2 - S3 - « 4 ) 

r(4 — D + Etl S *) T ( D -3- 3l -32- S 3 ) 

X T(3D/2 - 4 - £ti ^)r(3 + si + .S3 + s 4 - D/2) ' 
Closing the integration contours in a way that the convergence is guaranteed we get 

f -v 1 

\AlTIJL 



(23) 



Ti234(p 2 ]ml,ml,ml,ml) = f-r~^~2 ) <! > > B * + B ) ( 24 ) 



where 



B = {-X 4 ) 2V 1 p/^ + \ -^(1, 1 ~ V\ 1 + V, 1 - V\ Xi,X 2 ) 

[r 2 (^)r(l - 2v)F 4 (l - u, 1 - 2v\ 1 - v, 1 - v;r, t) 
+t u T(v)T{-v)F A {l, \-v\l-v,l + v;r,t) 
+r"T(v)T{-v)F A (l, l-v;l + v,l-v;r,t) 

^rfUl V) l + u;l + u,l + u;r, t) 



(25) 



with Xi = m 2 /p 2 with i — 1, 2, 3, 4, r — m\/m\ and t = m 2 /m 2 . Between the brackets we recognize up 
to some factors the result for the vacuum diagram with three massive lines || . Furthermore the coefficients 
c,; are 



Cl 



T 2 (u)T{u- l)r(l -2v) 



C2 = (-^2) 
C3 = (-^3) 
C4 = {-Xi) 



r(3i/-i) 

,r( i/)r(-i/)r(i/-i)r(2-»/) 
r(2v- 1) 

r(2i/- 1) 
„i»r(-i/)r(i-,/) 



r(2i/) 

c 5 = (-^n-^rr 2 ^), 
c 6 = (- X2 y(-x 3 yT 2 (-is). 
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The other terms are quartic series 





OO 

- E- 


OO 

■E 




ii — n 


J 4 — U 


B 2 


OO 

- E- 


OO 

■E 




3 1-0 


J4=0 


B, 


OO 

= E- 


OO 

•E 




31=0 




B A 


OO 

= E- 


OO 

■E 




31=0 


.7*4 =0 


B 5 


OO 

- E- 


OO 

■E 




ji—0 


J4=0 




OO 

- E- 


OO 

■E 




J i-0 


34=0 



1)31 (1 - - 2^3+34(2 - 



31+32+33+34 



(2 - ^ v )h+h+03+i 



(1 " ^32 (1 - ^33 (1 - ^34(2 - u) jl+h+ji (2 - 2v) h+js+ji 

l)3l(l _ V )j3+jA^ — 2^)j3+34(2 - 2^)^+^2+^3+^4(2 - V)j 1+ j 2+ j 3+ j i 



(1 + i/) j2 (l - i/) js (l - v) h {2 - v) h+h+h {2 - 2v) Jl+n+H 

1)31(1)33+34(1 ~ V )j3+3i(^ ~ +32+33+34 (2 — ^)j 1 +j 2 +j 3 +j 4 ^ 



(1 - v) h {l + i/) js (l - v) h {2 - v) 3l + n+H {2) n+n+H 

1)31(1)33+34(1 _ v )h+jiQ ~ 2^)31+32+33+34(2 - +32+33+34 



(1 - ^32(1 - V )js{ 1 + V )h ( 2 - ^)jl +33+34 ( 2 )3l+i3+34 
1)31(1)33+34(1 _ V )h+ji (2)jl +32 +33 +34 ( 2 ~~ ^)j 1+32 +33 +34 „ 



(1 + v)h (1 - u )h (1 + (2 - v)h+h+h (2)31+33+34 

l)3l(l)33+34(l ~ V )j3+ji (2)jl+J2+33+34(2 ~ V )jl+32+j3+ji 



(1 + ^32 (1 + ")j 3 (1 - ^34 (2 - ^31+33+34 (2)jl+j3+34 



where 



^,31 „32 „J3„J4 
1 2 3 4 

jib's bjtyj 



The large p 2 result given by (|24|) is valid for 

\p 2 \ > (rri2 + TO3 + mi) 2 and mi > TO3 + 7714 . 



(26) 



(27) 



(28) 



One may wonder what the relation is between the large p 2 expansion of (^4|) and that given in || . In 
the latter approach the various terms in the p 2 expansion are obtained from the expansion of subgraphs. 
The subgraphs are obtained by distributing the momentum p over the propagators in all possible ways. In 
the case of the diagram T1234 one has the following subgraphs: the diagram itself, the four diagrams where 
one internal line is removed, the two diagrams where two internal lines are removed and one diagram where 
three internal lines are removed. So one expects eight subgraphs, in fact seven because the contribution 
from the one where the line corresponding to propagator 1 in our convention is removed is zero. 

Following the analysis of || one can easily find the first term of each of the contributing series. For the 
subgraph representing the whole diagram the first term in the series should be the massless diagram. This 
series then corresponds to the first term in j2^), i.e. ~B\. The series which originates from the subgraph 
where two lines have been removed, e.g. 2 and 3, starts with the product of two massive tadpoles. They 
contribute a factor (m^rn 2 ^ 1 ' which can be identified with the sixth term in (pi]). The remaining subgraphs 
are obtained by removing one internal line, e.g. line 3. This yields a series starting with a massive tadpole 
proportional to {m\Y . This is the third term in (p^). The last term in ( p4| ) corresponds to the diagram 
where three internal lines are removed whose reduced graph is the vacuum diagram. Thus the seven series 
in (E4J) can be related directly to the seven subgraphs which are required for the method of 0. 
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The large p 2 expansion satisfies a system of four partial differential equations given below 

(T 1 {l + D 1 )x^ 1 -T 2 (l + D 1 ))f = Q, (29) 

((d + D 2 ){1 + D 2 )x 2 1 -T 2 )/ = 0, (30) 

(Ti(d + As)(l + AO** 1 - r 2(d + £» 3 + D 4 )(e + D 3 + Di)) f = , (31) 

(Ti(d + D 4 )(l + AK* - T 2 (d + D 3 + Di)(e + D 3 + D 4 )) / = , (32) 



where 



Ti = (a+ X) A)(c+ D *)> ( 33 ) 



i=l,3,4 i=l,3,4 
4 4 



r 2 = (a + ^A)(6 + ^A), (34) 

i=l i=l 

a = (2 - 2i/), & = (2 - 3i/), c = (2 - 1/), d = (1 - 2/), e = (1 - 2i/) and A = x l d/dx l . 

The expressions for small and large p 2 are given in an arbitrary number of dimensions and for arbitrary 
masses. One may subsequently derive special cases setting masses equal to zero. Taking the small p 2 
expansion, 7713 and TO4 cannot be taken zero simultaneously. When enough masses are zero and others are 
equal it is not so difficult to recognize known hypergeometric functions and to perform an expansion in 5. 
We give the following examples. 

We will use the following abbreviations 



^=7 + ln(— 5 , L p = 7 + ln( z £ 5 ) . (35) 



In the case mi = m 2 = TO3 = the three-particle cut contribution ( |l7| ) reads 

T£l 4 (p 2 ; 0, 0, 0, m 2 ) = Y{1-8)T{8-\)Y{28) 3 F 2 (1, «5, 28; 2, 2-5; x) (36) 
and the two-particle cut contribution is given by 

/ -p 2 \ ~ 5 ( m 2 \ ~ & r 2 (l- S)T{5)T{5 - 1) 
\4tth 2 ) \4tth 2 ) T(2-28) 



Bo(p 2 ;0,0)S (0;0,m 2 ) = -( 7 ^) (^) - ^ ^ - ^ . (37) 



Expanding the two results from above in 5 we obtain 

Ti 2 34(p 2 ;0,0,0,m 2 ) = -L + l||-i m -lrx(-^ 



-H - | C ( 2 ) + L 2 n + (-5 + 2 In (-*)) L m + 1 In 2 (-a;) 

-3 In (-s) - 2^—^ In (1 - x) - _-tlLi a ( x ) , (38) 



where x = p /m , in agreement with 
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In the case m\ = mi = and 7713 = TO4 = to the contribution from the three-particle cut ( |l7| ) becomes 

T (3] , 2^ o r(5-i)r(25)r(2+5) / m 2 \- 25 

r^(p) = 2 (^-j) 3*2(1, 5, 25; 2-5, 3/2+5; z/4) 

~i + ^G +im ) -^- L ™-^-^(2)+r + 0(5), (39) 



where 



S71 



T 2 (n) 



- (l + n)r(2n + 2)" 



-|( 3 F 2 (1, 1, 1; 2, 5/2; z/4) - i 3 F 2 (l, 1, 2; 3, 5/2; z/4)) . (40) 

Knowing the analytic expression for the Bq functions which occur in the two-particle cut contribution and 
expanding in S we obtain 

T 1234 (p 2 ;0,0,to 2 ,to 2 ) = 1^ + 1(1-^+^-^(2) 

+^-5L p + 31n(-f)-Iln 2 (-f 

m 2 , . . 1 / 2m 2 \ , 

+ 3 ^2-( r i - ^)ln(n) + - (l + ) ln 2 (n) , (41 ) 



where T\ and are the roots of the equation 

2 

m 2 r + — = 2m 2 - p 2 , (42) 
r 

in agreement with ||. 

When m 2 = TO3 = and mi = TO4 = to cq. (|l^) gives 



^) ^2_ 2 nnm 2,_ T(5)r(25)r(l-5) ( m, 

(l-5)(l-25) V 47r M 2 



2 x -2<5 



T i w 4 (^;m^0 ) 0,m^) = v ; ( ) 2^(25, 5; 2 - 5; x) . (43) 



To this the contribution Bq(p 2 ; m 2 , 0)_Bo(to 2 ; 0, to 2 ) should be added 

T^ 2 ; to 2 , 0, 0, to 2 ) = ^g^L ~" 2Fl(1 , 2 _ ff; ,) . (44) 

Expanding in <5 we get 

' 1 - x 



„2. m 2 



1 


1 


2P ^ 


" 25 


19 


1 


+ Y 


~ 2* 



X 

- 2 



ln(l - x) 
1 - .T 



C(2) + L m - 5L m - 2L m In (1 - x 
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+5 



In (1 - x) - 



In 2 (1 - x) 



-Li 2 (x). 



(45) 



In principle there is one class of cases which has not yet been covered by the above multiple series. 
That is the case where 7713 = 7724 = 0. For completeness we give the small p 2 case as a linear combination 
of an F2 and an F4 functions || [l(], |Tl[| . For the derivation we used the Mellin-Barnes representation 
method. 



Ti234(p 2 ; mf, mj, 0,0) 



-26 



r(i - s)r(s - i)r(2<5 - 1) 



( m l 

\4?T^ 2 

(1 ± ^Tj)- 2{1+5) F 2 (l + 6, 1,3/2 - 5; 1 + 5,3- 25; z ly z 2 ) 
-x^F^l, 2 - 5; 2 - 5, 2 - 5; x, y)] (46) 

where x — m\jm\, y = p 2 /mf, Z\ = x/(l±^/y) 2 and z 2 = ±4^/(1 ± y/y) 2 . 

As it was mentioned in the introduction the case T11234 with m[ = mi can be obtained from T 12 34 by 
differentiating with respect to m\. For instance, differentiating eq. ( |l9| ) leads to 

/ 2 2 2 2 2 2\ ^ (3) / 2 2 2 2 2\ 

Tii234(j> ;m 1 ,m 1 ,TO 2 ,m 3 ,m 4 ) = — ^T^{p im^^m^m^) 



dm 2 

+B Q (p 2 ;vi =2,i> 2 = hm^m^Bo^^m^ml 
+Bo{p 2 \m 2 l ,m 2 2 )-^^B (m\\ml,m 2 4 ) . 



(47) 



In this result the propagators in the B functions have the usual powers = 1 except when indicated 
v\ = 2. One can now even take the limit mi = and then expand the result in S. In section 4 more details 
on the products of Bq functions and their expansion will be given. 

From eq. (^) one can obtain special cases. Taking mi = and one other mass vanishing double series 
are obtained. In some cases one may obtain known hypergeometric functions, e.g. for mi = m 2 = 0, 
TO3 = m4 = to, which leads to 



11234 



(p 2 ; 0,0,0, m 2 , m 2 ) = 2 



47T/J 2 
2 



-26 



r(s - i)r(i + 2<T)r(3 + s) 

r(5 + 25) 



x 3 F 2 {l,S, 26+ 1; 2 -5, 6 + 5/2; 



4to 2 ' 



4?t^ 



-a 



/ 2 
/ — p 



\47r/x 2 



-6 



r(i + 5)rfflr 2 (i-5) 
r(i - 25) 



i 



i 



6m 2 (l - 25) p 2 5 



(48) 



3 Analytic approaches and elliptic integrals 

In this section we inspect the imaginary parts of the London transport diagram Ti 2 3 and of Xi 2 34. It turns 
out that they can be calculated in four dimensions in terms of complete elliptic integrals. These are well 
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known functions and thus the results arc of analytic interest. Furthermore fast and precise algorithms for 
the calculation of the elliptic integrals are available. Therefore the results provide also an efficient way to 
calculate the imaginary parts numerically. 



3.1 Imaginary part of the London transport diagram 

As can be seen from (^) the imaginary part of T123 is convergent in four dimensions and reads with a 
factorization of the Kallen functions 

1 . ™ , ■)> 7T f dt 



Im{T 123 {p 2 )) = -AT 123 (p 2 ) = ~j jy/(t - xi){t - x 2 ){x 3 - t)( Xi - t) , (49) 

with 

X\ = (mi - m 2 ) 2 ; x 2 — (mi + m 2 ) 2 ; x 3 = (p - m 3 ) 2 ; X4 = (p + m 3 ) 2 ] 
X\ < x 2 < x 3 < X4; p = \/\p 2 \ > mi + m 2 + m 3 . 



The integration limits are zeros of the square roots, and thus (49) leads to complete elliptic integrals, 
defined by 

1 

f dt 7T „ , 1 1 0, , 

K(x) = / m = ==^- 2 Fi(~-,-;l;x 2 ), (50) 





1 



y/(l-t 2 )(l~X 2 t 2 ) 



dt(l-X 2 t 2 ) TT 11 



E W = / m K == ±= = -2F 1 (-,~;l;x 2 ), (51) 





1 



y/(l -t 2 )(l ~X 2 t 2 ) 



2 v 2 2 



ll(c,x) = I = = ~-Fi(~;l, \;l;c,x 2 ), (52) 

' (l-ct 2 )^(l-t 2 )(l-x 2 t 2 ) 2 u 2' '2' ; ' v ' 



with the Gauss hypergeometric function 2 F\ and the Appell function F\ || |T^, pT| . Reduction of ( fl9| ) 
to the Legendre normal form of the elliptic integrals [[12| |l3j by decomposition into partial fractions and 
partial integration yields after some algebra 

7m (Ti23(p 2 ;m 2 ,ml,m 2 )) = 



-)2 



' Amim2[(p + m 3 ) — m 3 p + mim 2 \. K 



q++ VV 



ml+ml + ml+p 2 / q+-q- 

2 VV 

8mim 2 [(m 2 + mjjfp 2 + m§) — 2m 2 m| — 2m 2 ; p 2 ] ^ g_ + jq + _q_ + 



y/q++q — V<7 — V 

8mim 2 (p 2 - m§) 2 ^ (mi - m 2 ) 2 jq + -q- + 



^/q++q — V(mi+m 2 ) 2 g — \ q++q- 

xO(p 2 - (mi + m 2 + m 3 ) 2 ) , (53) 
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with variables q±± corresponding to the physical and unphysical thresholds 

q±± ■■= (p ± m 3 ) 2 - (mi ± m 2 ) 2 ■ (54) 

This result is valid in all parameter regions. In special cases it leads to simpler formulae. For equal 
masses one gets 

Im(Ti23(p 2 ;m 2 ,m 2 ,m 2 )) = 7 J (p - m)(p + 3m) 

x{ - Am 2 p K{k) + (P~H(p 2 + 3m 2 ) ^ j Q( ^ 2 _ ^ ^ 

.,, 2 (p + mf(p-3m) 

with k := - , (56) 

(p — m)°[p + 6m) 

involving only complete elliptic integrals of the first and second kind, i.e. 2-F1 Gauss hypergeometric 
functions. If at least one mass is zero, Im(T\23) reduces to logarithms. The most complicated case leads 
to 

r ,rr i 2 n 2 2nn n f p 2 + m 2 +m 2 

Im(Ti23(p ;0,m 2 ,m 3 )) = --^ j y/q+q- 



2 2 2 , 2 , 2\\ 1 ,P + m 3 - ™2 - 



+ (2m^ 2 -m^(m|+p 2 ))log(- 



2m 3 p 

+m\(p 2 - ml) log( — — 2m 2 m 3 p — ) 

(jr — m 3 ) z — m 2 (?7J 3 + p ) + (p — m^)Jq+q- 

xO(p 2 - (m 2 + m 3 ) 2 ) , (57) 
with g± := (p ± m 3 ) 2 — m\. 

3.2 The imaginary part of T1234 

The two-particle cut contribution to the discontinuity of T1234 was given in ([l3]) , 

A r i2L(p 2 ; TO ?> m 2> 7^3,7714) = AB (p 2 ;ml,ml)B a (ml+ie;ml,ml) . (58) 

As a product of a one-loop self-energy integral and a one-loop self-energy discontinuity it is composed of 
elementary functions and gets a real part for 

(m 3 + m^) 2 < m\ and (mi + TO2) 2 < p 2 . (59) 

The three particle cut contribution can be calculated in a fashion very similar to the case of the London 
transport diagram. Only one more (complex conjugated) propagator l/(t — m 2 — ie) has to be added in 
(|H). The calculation yields 

at-t(3) / 2 2 2 2 2\ 
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P 2 I V 1++ \\ Q++1— J \ V 



n 



\3— V <?++<?-- 
8m 3 TO4(7J 2 - m 2 .) 2 ^ (m 3 - 777 4 ) 2 q_ + /<7 + _<7_ + 



m iy / 1++1 — \{m 3 + m 4 ) 2 q — ' \ q++q- 

| 87773777 4 A(p 2 ,777 2 ,777^) ^ / 7Ttf — (TO3 - TO 4 ) 2 g-+ ^ 

m iy / q++<i — - ( TO 3 + »»4) 2 g — ' V '?++'? — 

x9(j) 2 - (7773 +7773 +TO4) 2 ) , (60) 



with 



<7±± := (P± "72) 2 - (m 3 ± 777 4 ) 2 . (61) 
In the case (^9|) the characteristic c of the last II- function in ( |60| ) is greater than 1, 

w 2 - (7773 ~ w 4 ) 2 q-+ 
i\ — (7773 + m 4 ) 2 q- 



c= 2 LM 4 ( 2 F L >1, (62) 



which requires an analytic continuation of that function. A comprehensive discussion of the analytic 
properties of the elliptic integrals can be found in [ p^[ . The m 2 + ie prescription in ( |60| ) ensures that 
AT-^234 gets the correct real part, given through 

Im (II(c — ie, k)) = — (II(c — ie, n) — H(c + ie, k)) 
2i 



= _£/ 1 _. (63) 

2 V (c- 1)(c-k 2 ) V ' 

This contribution cancels the real part of the two-particle cut AT 12 3 4 - Consequently AT1234 is always 
purely imaginary. 

Numerical checks show the agreement of the results of (|53|) for Im(Ti23j and of 



Im(T 1234 ) = i (aT x ( 2 2 34 + AT x ( 2 3 y (64) 



with previously published tables 



4 One-dimensional integral representations 

4.1 A general approach to two-loop integrals containing a self-energy subloop 

An alternative method to the series expansion of the two-loop scalar diagrams consists in the derivation 
of one-dimensional integral representations. These are built up from one-loop self-energy functions Bq 
coming from the self-energy subloop and the remaining one-loop integral. They can be derived by using a 
dispersion representation of the B function. 
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A two-loop diagram with only three-vertices 



PN 




where k is the momentum flowing through the self-energy insertion, can in a first step be reduced to 
simpler diagrams by a decomposition into partial fractions 



k 2 m 2 N+2 k 2 m ?V+2 m 7V+3 m JV+2 k 2 m %+3 

This yields for the diagram 



m N+2 m N+3 



Ti...N + z{pi,m 2 ) = — — g — ( Ti...jv+2(Pi;mi, . . . ,m% +1 ,m 2 N+2 ) 



-Ti... N +2{Pi;m\, . . . ,m 2 N+1 ,m 2 N+3 )j 



The difference has to be replaced by a derivative if rn 2 N+2 ~ m ?v+3- 

Insertion of the dispersion representation for the self-energy subloop leads to 

T 1 ... N +2(Pi;m 2 1 , . . . ,m 2 N+1 ,m 2 N+2 ) 
= ( B o(k 2 ;m 2 N ,m 2 N+1 ) {k+p l )2 _ mV -- 
1 1 



(k +px + . . . -\-pN-i) 2 - k 2 -m% +2 

oo 

!/".._. o « . / -1 



— / dsAB (s;m^ra 2 w+1 )/ — 



so 



with So = (mjv + ?^7v+i) 2 - After a further decomposition into partial fractions, 

-1 11/1 1 



(65) 



X fc 2 - (fc +Pi) 2 - m 2 " ' (fc +pi + . . . +Pat-i) 2 - m 2 v _ 1 ^ *~ 66 ^ 



(67) 
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the fc-integrations and one of the s-integrations can be performed and yield 

Tx...N+2(pi', mj,..., m 2 N+1 ,m 2 N+2 ) 
= B^{m 2 N+2 ]m 2 N ,m 2 N+1 )T { - 1) {p l -,m\, . . . , m^_ lt m 2 N+2 ) 



1 [°° AB Q (s,m 2 N ,m 2 N+1 ) (1) 2 2 

— / ds — 3 T K >(jp l ;m 1 ,...,m N _ 1 ,s). (68) 

s s m N+2 



2ni 



denotes a one-loop N-point function in which s enters in the remaining one-dimensional integration 
as a mass variable. 

A diagram with two four- vertices leads to a result which is similar to the remaining integration in ( |68|) , 



jn N -i 

Pn-i 




m N+ i 77i at 



OG 

Ti... N+1 (pi;m 2 ) = — dsAB a (s;m 2 N ,m 2 N+1 ) 
2m J 



k 2 - s (k +pi) 2 - mf" ' (k + pi + . . . +pn-i) 2 - ■ ,,J 



N-l 



~~2~K~i J dsAjB o(s; 777^,m^ +1 )T (1) (p l ;777i,...,m^_ 1 ,s). (69) 



4.2 Examples 

An application of (|69j) to the London transport diagram leads to 



OG 

Ti23(p 2 ;m 2 ,ml,ml) = j ds AB Q (s;ml,ml) B o(p 2 ; s,m\) , (70) 

(m 2 +m 3 ) 2 

a result which would also follow from (||). In that case a suitable subtraction || is 

T 12 3n(p 2 ; m\, 777,3, m l) = Ti23(p 2 ;m 2 1 ,m 2 ,,ml) -Ti23(p 2 ;ml,0,ml) (71) 

-Ti23(p 2 50,7772,7773) +Ti 23 (p 2 ; 0,0,7773) . 

For X1234 one obtains from ( |68| ) 

T1234 (p 2 ; m 2 , , 7773 , m 2 ) 
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1 

2Tri 



AB (s;ml,ml) , 2 2 2 2 2 . 

ds 2— — (£ (p ;m 1 ,m 2 ) - Bob ; s ) TO 2)j 



(m3+m 4 ) 2 
,2. 2 2 



B (mi;m 3 ,m 4 )So(p ;m 1 ,m 2 



1 



A-B (s;m§,ml) 2 2 
2— B (p ;s,m 2 ) 



(72) 



(73) 



(rri3+m 4 ) 2 



The representations ( JTOj ) and ( [73] ) with the subtractions ( J7l] ) and (|2(j) provide efficient ways to calculate 
T\23N and Ti 2 34at in all parameter regions. The results agree numerically with those published in jj). 
One may also consider vertex functions, for example 



P3 



-Pi - P2 




rrl 222 2 2\ o/2 2 2\,^(/222 2 2 2\ 

T{p 1 ,p 2 ,P3; m iT--, m 5) = B a{m 1 ;rn 4 ,m 5 )Co{p 1 ,P2,P3;m 1 ,m 2 ,rn 3 ) 



1 

2tH 



AB {s;ml,m 2 5 ) 2 2 2 2 2 

ds 2 Co(p x ,i> 2 ,i>3; s >' Tl 2j"' l 3). (74) 



In that case the counterterm of the self-energy subloop can be subtracted for a numerical evaluation. 

In the cases of Ti 23 jv and Ti 2 34jy equivalent one-dimensional integral representations can be obtained 
with a direct integration in the momentum space or using the dispersion representation. 

The case Tn 2 34 for m x = gets the representation 



7ii234b 2 ;0, 0, ml, ml, ml) = B (p 2 ;vi =2,l> 2 = 1; 0, m\)B {Q; m\, ml) 

00 

1 f , AB (s; ml, mf) 



2ni 



As 



[B (p 2 ;s,m 2 2 ) - B (p 2 ;0,ml)] . 



(75) 



The integral is convergent and the result in an arbitrary number of dimensions for the product BqBq which 
we denote by Z is given by 



-1 / 9 \ - ^ / ? 
1 / m 2 \ / m\ 



% = 



-6 , 2 x -1 

TO4 



1 - 



l-<5' 
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^2 

One can perform the expansion of Z in <5 and obtains 



x H7T + ,? 2 2^1(1 + 2; 2 - ftp 2 /™!)- (76) 



+a (L_ 4 + b - 2) + \ (L m4 + b) 2 -2 [L mi + b) - ^| [_- 4 -► _ m3 ] }> , (77) 

Z TOyl 



'4 



where a; = p /to 2 ,, a = 1 + (1/cc — 1) ln(l — x), 6 = L r „ 2 + 21n(l — x) and 

X = f — ln(l -x)-( —) ln 2 (l - x) + f — ') Li 2 (x) . (7.M 



In the above representation masses can be set to zero, except the case 1713 = 7774 = 0. However, for this 
case an explicit result can be derived 

T 11234 (p,0,0,TO,0,0) = _(^ — j _____ 

x 2 Fi(l + 25, 2 + 5;2-S;x) (79) 
with a; = p 2 /m 2 . Expanding the result around D = 4 we get 

T 11234 (p 2 ;0,0,TO 2 ,0,0) = ^ 2T b^){^-K C + rf J) +9 + 3C(2) 

+2d 2 -6d + 4c(d-^J+4 (^r^j ln2 (! - - T ) 

' " n in(l-.x) + 2^^±l)Li 2 (x)| , (80) 

where c = 1 + — 1) ln(l — a;) and d = _ m + 21n(l — a;). 

For the other cases where some masses are zero the integral can also be calculated analytically. A 
simple case is 

P 2 T 11234 (p 2 ; 0, 0, 0, m 2 , 0) = ~ + i {-1 + L p + L m } - | - i (_ p + L m f 

+L P + L m - -a; In (-x) + Li 2 (x) 

+ i(x-i)ln(l-_) (81) 

where a; = p 2 /ttt, 2 , in agreement with M. 
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5 Numerical comparisons and conclusions 



In this section numerical comparisons are made between the various results of this paper and the literature. 
Moreover some concluding remarks and an outlook are presented. Various numerical evaluations for 
small p 2 are presented in Table 1. In the first column are values from the two-dimensional numerical 
integral representation in the second results from the series obtained from the hypergeometric functions, 



(appendix), whereas column three uses three terms of the Taylor expansion of 15 and finally the last 
column is calculated from the one-dimensional integral of sect. 4. The comparison gives confidence in the 
various calculational methods. 

In conclusion the analytic results in terms of generalized hypergeometric functions offer the possibility 
to use well-known mathematical techniques for analytic continuation, partial differential equations, contour 
integral representations, etc. Moreover the result is derived for arbitrary masses from which one can get 
many special cases. In other words the several formulas for the cases with vanishing masses are unified 
in one result in D dimensions. When one is only interested in the imaginary parts an alternative analytic 
result in four dimensions is obtained in terms of complete elliptic integrals. Since their properties are 
well-known they are easily accessible for numerical evaluations. Finally we have derived a one-dimensional 
integral representation for all two-loop diagrams containing a self-energy insertion. For the two-loop self- 
energy diagrams treated in this paper the integrand is composed of elementary functions only and the 
representation is valid for all values of p 2 . The main application of these integrals in this paper is for 
numerical evaluations, giving a good alternative to the existing two-dimensional integrals. In order to 
complete the treatment of massive two-loop self-energy diagrams one has to make the similar study of 
the master diagram. As to the analytic approaches the prospects are not so good since it is not of the 
self-energy insertion type. In practice this means for instance that the massless case with arbitrary powers 
of the propagators is not available in the literature. This prevents an application of the Mellin-Barnes 
representation. On the other hand symmetry properties, special cases and imaginary parts have been 
studied in jl6| . The two-dimensional integral representation was derived in particular for this diagram |jj . 
Nevertheless further studies will be needed. Once adequate techniques for the self-energy diagrams are 
available one could envisage practical applications for physics predictions. For the electroweak theory the 
obvious application is to the gauge boson self-energies which play a role in the My/ — Mz mass relation 
and details of the Z line shape. For QED the two-loop vacuum polarization has been known for a long 
time @, but the electron two-loop self-energy was never fully calculated. 
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Appendix 

Here we give the result for the small p 2 expansion in 5. We introduce the dimensionless variables x, y, u 
and v defined by 

X =^,y=!4, (82) 



20 



P 2 


A 


B 


C 


D 


2.66667 


-8.454752 


-8.45044 


-8.445577 


-8.45038 


1.77778 


-8.286016 


-8.28753 


-8.286116 


-8.28747 


1.18519 


-8.186514 


-8.18486 


-8.184417 


-8.18481 


0.790124 


-8.116017 


-8.11883 


-8.118665 


-8.11878 


0.526749 


-8.077571 


-8.07583 


-8.075740 


-8.07577 


0.351166 


-8.053491 


-8.04759 


-8.047528 


-8.04754 


0.234111 


-8.033367 


-8.02896 


-8.028900 


-8.02890 


0.156074 


-8.013727 


-8.01662 


-8.016561 


-8.01656 


0.104049 


-8.018667 


-8.00843 


-8.008371 


-8.00837 


0.069366 


-7.999787 


-8.00298 


-8.002926 


-8.00292 


0.046244 


-7.997535 


-7.99936 


-7.999303 


-7.99930 


0.030829 


-7.985865 


-7.99695 


-7.996891 


-7.99689 


0.020553 


-8.006871 


-7.99534 


-7.995285 


-7.99528 


0.013702 


-7.900682 


-7.99427 


-7.994214 


-7.99421 


0.009135 


-7.957410 


-7.99356 


-7.993501 


-7.99350 


0.006090 


-7.793835 


-7.99308 


-7.993026 


-7.99302 



Table 1: Real part of Ti234jv(p 2 ; m\ , m%, m§, mf) for small values of p 2 . The masses are mi = 1, m,2 = 3, 
TO3 = 5 and tua — 7. A: Result from VEGAS integration . B: Result of eqn. (B6J). C: Approximation 
using 3 terms of Taylor expansion. D: One-dimensional integral representation (h 



and 



mf mi 

u = ^,v = —1. (83) 
mi m 



'4 



We also need the roots r\ and of the equation 



2 



r 



and the roots r% and of the equation 



m\r H = ml + mi - m\ . (85) 



where 



Ti234aKp 2 ; ml,ml,ml,ml) = -Al + A2 + CI + C2 + C3 

+C4 + C5 + C6 + CI (86) 



Al = 8 -2C{2) + Ul+ V -^) hSuj) 
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+ (-2 (l + y—^j (In(n) - Info))) Info) 

1 — x , , , 1/1 — v + a;\ T .., 
+ --ln(l-x) + -^^^jL l2 (x) 

+ ^ ! ^{ 41n(ri) " 41n(r2) 
+ Li 2 f^lMVLi 2 ^ 1 -^ 



ri - r 2 y \ r 2 - ri 

-Li 2 ( -1^-) + Li 2 f — ^ - Li 2 (l - n) + Li 2 (l - r 2 ) 

1*1 — r 2 / \ r 2 — 1*1 J 

+Li2 ,ii(i^)u,ri-x,iv^' (i - r2: 



-r 2 / \ r 2 y \ -r 2 



A2 = |-2 + Info) - (Info) - In(r 2 )) 

2 + Y^r Hv) y 1 ^ 1 (ln(r3) ~ ln(r4)) 
+ \{~ 2 + tit Hy) ~ t 1 -^ 1 (ln(ri) ~ ln(r2)) 



+ l{-2 + I^ln(,)-i^(ln(r 3 )-lnM) 

1 , 2 m s 1 , 2 m ? 1, m? , to? 1 , m? , 

+ a lo § -4 + a l0 § ~4 + 7 lo S ~4 lo § "4 + 7 lo § ~4 lo § 

8 T114 8 m2 4 m| TO2 4 TO4 

+ ^{c(2) + 8+iln 2 ( y ) 
V ~ 1 Info) + r .±—^ (2 Info ) - 2 ln(r 2 ) 



2 


+ log 


9 

m| 
> 2 


+ log 


TO 2 






TO 2 




to| 


2 




9 








+ log 




+ log 


m| 


m| 




; to! 




m| 






2 






to| 


+ log 












TO| 







a; a; 



ri-r\) \ n - r 2 / V r i _r 2/ V r 2 - ri 

+Li2 , fiil^ _ Lis (aSl^ _ Li! C ^ + Li / 1 1 ! 



ri-r 2 / \ r 2 -ri 1 \ri-r 2 / \r 2 -ri 



+ i|c(2) + 8+iln 2 («) 
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-2- — -In(«) + - — — ( 2In(r 3 ) - 21n(r 4 ) 

+ lnfl^Vnf^ 1 -^Vlnfl^Vn 
\ r 4 - r 3 J \ r 3 -r 4 J \r 3 -nj 

+Li 2 (!±L^1) _ Li2 f ) _ Li2 (_LJl) + Li2 f lz£» ] u . 



r 4 (l - r 3 ) 



OO 



I _ . (m + n)!(m + n- l)!(m + n + fc + Z- l)\(m + n + k + l)\ 
to!(to + l)\n\{n - l)!/!(2(m + n + k) + I + 1)! 

fe';=o 

x [--0(m + n) — 2^>(m + n + fc + I) - ip(m + n + k + I + 1) 

-V(m + 2) + ^(n) + 2V>(2(m + n + fc) + 2 + /)] - z 3 ) l z?. (89) 



^ 2 = E 



(to + n)\(m + n + l)!(m + n + fc + / + l)!(m + n + fc + /)! 



m!(m+l)!n!(n + l)!/!(2(TO + n + fc)+Z + 3)! 
x [log(ml/ml) + ip(m + n + k + I + 1) + tp(m + 2) 

+V>(n + 2) - V(m + n + 2)] - z 3 ) l z?. (90) 

^ (m-l)!(m + fc + Z-l)!(m + fc + Z)! fc , 
C3 ^ - fe ^ (m + m{ 2(m + k) + l + l) (1 - * 3 > *< (91) 

m — 1 

C4 " fe ^o l!(2(n + *) + / + !)! 'i^ 1 "^ 

rl=l 

[-1 + 2^(2(n + fc) + I + 2) - 2^(n + k + I) - ip(n + k + I + 1)] . (92) 

^ (k + I - iy.(k + l)\ kM .,. 1 
= E / , (2fc + ^(l-^)'[-l-^(fc + 

+ Z + 1) + 2-0(2fc + I + 2)] . (93) 

06 = E Sil 'i t- 1 - 2 ^( fc ) - ^( fc + !) + 2 ^ 2fc + 2 )i ( 94 ) 
ct = E frrf (1 _ Z3)l [_1 ~ 2?K0 _ ^ + 1} + 2Vj(/ + 2)] (95) 
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